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Abstract 

We present Keplerian-type parametrization for the solution of third post-Newtonian (3PN) ac- 
curate equations of motion for two non-spinning compact objects moving in an eccentric orbit. The 
orbital elements of the parametrization are explicitly given in terms of the 3PN accurate conserved 
orbital energy and angular momentum in both Arnowitt, Deser, and Misner-type and harmonic 
coordinates. Our representation will be required to construct post-Newtonian accurate 'ready to 
use' search templates for the detection of gravitational waves from compact binaries in inspiralling 
eccentric orbits. Due to the presence of certain 3PN accurate gauge invariant orbital elements, 
the parametrization should be useful to analyze the compatibility of general relativistic numerical 
simulations involving compact binaries with the corresponding post-Newtonian descriptions. If 
required, the present parametrization will also be needed to compute post-Newtonian corrections 
to the currently employed 'timing formula' for the radio observations of relativistic binary pulsars. 

PACS numbers: 04.30.Db, 04.25.Nx, 04.80.Nn, 95.55.Ym 
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I. INTRODUCTION 



An accurate analytic description for the orbital dynamics of compact binaries of compara- 
ble masses is highly desirable to probe astrophysical scenarios involving strong gravitational 
fields. These include inspiralling of compact binaries, relevant as sources of gravitational 
radiation for the operational and proposed interferometric detectors and relativistic binary 
pulsars. 

The method of matched filtering, ideally suited to detect gravitational radiation emit- 
ted by inspiralling compact binaries, requires 'ready to use' search templates These 
templates usually consist of gravitational wave polarizations, h + (t) and h x (t), temporally 
evolving due to the orbital dynamics of the binary. In the inspiral regime, the orbital dy- 
namics is well described by the post-Newtonian (PN) approximation to general relativity, 
which allows one to express the equations of motion for a compact binary as corrections to 
Newtonian equations of motion in powers of (v/c) 2 ~ GMj (c 2 R), where v, M, and R are the 
velocity, total mass and relative separation of the binary. Currently, widely employed search 
templates are for compact objects moving in quasi-circular orbits, where h + (t) and h x (t) are 
known, both in amplitude and phase evolution, to second post-Newtonian (2PN) order Q]. 
However, it is highly desirable to know a priori the temporal phase evolution of h + (t) and 
h x (t) at least up to third post-Newtonian (3PN) order beyond the (Newtonian) quadrupole 
contribution which in turn requires the knowledge of the conservative binary dynam- 
ics to 3PN order, i.e, (v/c) 6 corrections to Newtonian equations of motion. Very recently, 
after tackling for years many conceptual and computational issues, the phase evolution of 
gravitational wave polarizations to 3PN order was achieved, [see jj] and references therein] . 
Also recently, a formalism to compute 'ready to use' search templates for compact binaries 
in inspiralling eccentric orbits was developed and explicitly implemented at 2.5PN order 
Its natural extension to higher PN orders, which should be relevant for detecting compact 
binaries in eccentric orbits, requires a parametric solution to the 3PN accurate equations of 
motion. 

To describe the late stages of binary inspiral and the subsequent merger, the emphasis is 
currently being placed on numerical relativity,which attempts to solve the associated set of 
full Einstein equations using supercomputers [6|. However, these general relativistic simula- 
tions usually do not incorporate inputs from PN dynamics and are incapable of describing 
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few binary orbits even in the inspiral regime. Recently, a 'post-Newtonian diagnostic tool', 
was introduced with the aim of extracting some physical information from the existing gen- 
eral relativistic simulations [?], This 'diagnostic tool' requires a solution of 3PN accurate 
equations of motion for compact binaries moving in non-circular orbits. The tool further in- 
troduced certain definitions for the eccentricity and semilatus rectum parameters in terms of 
the orbital angular velocity at the turning points of the orbit, making many of its estimates 
gauge dependent. It is desirable to characterize a non-circular orbit in post-Newtonian rela- 
tivity in terms of gauge invariant quantities and the parametrization we present here should 
be useful in that aspect. 

Finally, we note that the high precision radio- wave observations of binary pulsars employ 
an accurate relativistic 'timing formula' fl Q which requires an analytic solution to the 



relativistic equations of motion for a compact binary moving in an elliptical orbit . This 
'timing formula' is instrumental to test both the predictions of general relativity and the 
viability of alternate theories of gravity in strong field situations Q]. It should be noted 
that the long term measurements of the recently discovered relativistic double pulsar system 
J0737-3039 will require the inclusion of higher order PN effects in the timing formula jl^ . 

In this paper, we derive an analytic parametric solution to 3PN accurate conservative 
equations of motion for compact binaries moving in eccentric orbits. The orbital representa- 
tion is given both in Arnowitt, Deser, and Misner (ADM)- type and harmonic coordinates. 
We employ similar techniques which allowed Damour and Deruelle to obtain a remarkably 
simple parametrization for the solution of 1PN accurate equations of motion for compact 
binaries in eccentric orbits Further, our 3PN accurate orbital representation is struc- 
turally quite close to the generalized quasi-Keplerian representation obtained by Damour, 
Schafer, and Wex for the solution of 2PN accurate orbital motion of compact binaries in 



eccentric orbits in ADM gauge 



orders, a slowly precessing eccentric orbit can be characterized by certain gauge invariant 
quantities. This feature of the parametrization should make it very attractive to describe 
eccentric orbits in post-Newtonian relativity. We note that apart from some misprints, an 
incomplete representation in ADM- type coordinates was obtained earlier 

We have the following plan for the paper. In Sec.|nJ we first review Keplerian parametriza- 
tion associated with the Newtonian accurate orbital motion. This will be followed by brief 
descriptions about the quasi-Keplerian and the generalized quasi-Keplerian parametriza- 
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tions associated respectively with the 1PN and 2PN accurate orbital motion. Sec. IIHI deals 
with the derivation of 3PN accurate orbital representation for compact binaries moving in 
eccentric orbits in ADM-type coordinates. In Sec. II VI a similar parametrization, valid in 
harmonic gauge, is derived. Finally, in Sec. we summarize our results and discuss the 
merits and applications of the parametrization. 



II. THE KEPLERIAN PARAMETRIZATION AND ITS PN EXTENSIONS 

The Keplerian parametrization for Newtonian accurate orbital motion of a binary in 
eccentric orbit is heavily employed in celestial mechanics [l^. In polar coordinates and in 
the center-of-mass reference frame, the eccentric motion is parametrized in the following 

way 

R = a(l — ecosu), (la) 



2 arctan 



, 1/2 

e \ u 
tan — 



(lb) 



l-ej 2 

where R and <fi define the components of the relative separation vector R = i?(cos 0, sin 0, 0). 
The semi-major axis and the eccentricity of the orbit are denoted by a and e respectively; 
both are expressible in terms of the Newtonian conserved (orbital) energy and angular 
momentum. The auxiliary angles u and v are called eccentric and true anomaly respectively. 
The geometrical interpretation of these anomalies is presented in Fig. ^ The explicit time 
dependence is provided by the Kepler equation, which reads 

I = n{t — t ) = u — e sin u , (2) 

where I is the mean anomaly and n is referred to as the mean motion and is given by n = 
P being the orbital period. The quantities to and <fro are some initial time and orbital phase. 
The explicit functional dependence of the orbital elements a, e and n is given by 

GM 

a = (-2~B) ' (3a) 
e 2 = l + 2£/i 2 , (3b) 
(-2E) 3/2 , . 

where E is the Newtonian orbital energy per unit reduced mass \x = mim^/M , m\ and 
m 2 being the individual masses of the binary and M = mi + m 2 . The reduced angular 



momentum h is given by h = (jjj, where J is the reduced Newtonian orbital angular 
momentum. 

For 1PN accurate equations of motion, in harmonic coordinates, Damour and Deruelle 
found the following 'Keplerian like' parametrization which reads 



R 

I = n(t — to) 
^ (0 " 0o) 



a r (l — e r cosu) 
u — e< sin u , 

v = 2 arctan 



1 + 



1/2 



U 

tan — 
2 



(4a) 
(4b) 

(4c) 



Note that the three eccentricities e r ,et and e^, which differ from each other by PN cor- 
rections in terms of E and the finite mass ratio r] = fi/M, were introduced such that the 
parametrization looks 'Keplerian' even at 1PN order. The factor ^ gives the angle of ad- 
vance of the periastron per orbital revolution. Due to these features, in the literature, the 
above representation is usually referred to as the 'quasi-Keplerian' parametrization for 1PN 
accurate orbital motion. The parameters a r , e r , e^, n and e t are some 1PN 'semi-major' axis, 
'radial eccentricity', 'angular eccentricity', 'mean motion' and 'time eccentricity' respec- 
tively. These orbital elements depend on the 1PN accurate conserved energy and angular 
momentum; their explicit dependence was derived in 

The above orbital parametrization was extended, in ADM coordinates, to 2PN order by 
Damour, Schafer, and Wex Q. The 2PN accurate orbital pa—ation ha 8 the 
following form 



R 

I = n (t — t ) 
^ (0 - <h) 



a r (1 



e r cosu 



u — e t sin u 



9u\ 
c 4 ) 



[V — U) 



fat 
c 4 



smt> 



v + I — 7- sin 2v 



(5a) 
(5b) 

(5c) 



where v = 2 arctan 



1— e d 



1/2 



tan | 



The presence of the orbital functions g it , fat,9i<f> an d fa<f> indicates that the structure of 
the solution is more general than the 'quasi-Keplerian' one, prompting to coin the above 
parametrization as the 'generalized quasi-Keplerian' parametrization for 2PN accurate or- 
bital motion for compact binaries in eccentric orbits. The explicit expressions for the orbital 
elements <3>, along with g iu fat, g^ and fa^ in terms of the conserved 2PN ac- 
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curate energy and angular momentum as well as the finite mass ratio were derived, in ADM 
coordinates, in jis| . 

In the next two sections, we will derive a similar orbital parametrization for both ADM- 
type and harmonic gauges, which will analytically describe the solution to the 3PN accurate 
equations of motion for a compact binary moving in an 'eccentric' orbit. 



III. THE 3PN ACCURATE GENERALIZED QUASI-KEPLERIAN REPRESEN- 
TATION IN ADM-TYPE COORDINATES 



In the Hamiltonian formulation of general relativity, advocated by Arnowitt, Deser, and 
Misne r llal. an ordinary Hamiltonian for a compact binary is achievable only up to 2PN 
order (ijj, |20| . 

At the third post-Newtonian order, an higher order Hamiltonian which depends on par- 
ticle positions, conjugate momenta and their derivatives was derived by Jaranowski and 
Schafer j^]. Later, using an higher order contact transformation, the above Hamiltonian 



22] 



was transformed into an ordinary Hamiltonian by Damour, Jaranowski, and Schafer 
We refer to this ordinary Hamiltonian as the 3PN Hamiltonian in ADM-type coordinates. 
The two unknown coefficients, which initially appeared in the above Hamiltonian, were later 
fixed (23I. We display below the fully determined reduced ordinary 3PN Hamiltonian, in 
ADM-type coordinates and in the center-of-mass frame, compiled using 23[ 2^| as 



H(r, p) = W (r, p) + ^Wi(r, p) + ^H 2 (r, p) + -H 3 (r, p) 
where the Newtonian and post-Newtonian contributions are given by 



W (r,p) 
Hi(r,p) 
n 2 (r,p) 



n 3 (r,p) 



p_ 

2 
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-2r/ 2 (n-p) 2 p 2 -3r/ 2 (n-p) 4 
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(5 + 8r/)p 2 



1 1 

r + 2 



+3?7(n ■ p) 2 
1 



r 2 ^ 



2\4 



128 



(-5 + 35r/ - 70r/ 2 + 35r/ 3 ) (p 2 ) 



(6) 



(7a) 
(7b) 



(7c) 



1 

16 



2^ 



(-7 + 42r/ - 53r/ 2 - 5r/ 3 ) (p 2 ) 



(2 - 3r/)r / 2 (n • p) 2 (p 2 ) 2 + 3(1 - rj)r) 2 (n ■ p) 4 p 



4-2 



— 5?7 3 (n ■ p) 



2\2 



16 



(-27 + 136r/ + 109r/ 2 ) (p 2 ) 



12 
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I 192 
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+ ^[12+(872 - 63tt 2 ) 



(7d) 



where r = R/(GM), r=|r|,n = r/r and p = P /(/,', R and P are the relative separation 
vector and its conjugate momentum vector. 

The invariance of Ti under time translation and spatial rotations leads to the following 
conserved quantities: The 3PN reduced energy E — Ti and the reduced angular momentum 
J = r x p of the binary in the center of mass frame. The conservation of J particularly 
implies that the motion is restricted to a plane and we may introduce polar coordinates such 
that r = r(cos</>,sin0). The Hamiltonian equations of motion, which govern the relative 
motion read 

an 

dp 

dn 



n 



r 2 (j) 



r x 



dp 



(8a) 
(8b) 



where r = dr/dt,4> = d<j)/dt and t denotes the coordinate time scaled by GM. We intro- 
duce s = 1/r so that r 2 = (ds/dt) 2 /s 4 = s 2 / s 4 and obtain, using Eqs. (jBJ), 3PN accurate 
expressions for r 2 and d<j)/ds = 4>/s, in terms of E, h — | J|, rj and s, which are displayed in 
Appendix EI We note that the 3PN accurate expressions for f 2 and <fi are seventh degree 
polynomials in s. 

To obtain the 3PN accurate orbital parametrization, we proceed as follows. First, we 
concentrate on the radial motion and compute the two nonzero positive roots, having finite 
limits as - — * 0, of the 3PN accurate expression for r 2 . These two 3PN accurate roots, 
labeled s_ and s+, correspond to the turning points of the radial motion and hence to the 
periastron and the apastron of the post-Newtonian eccentric orbit. With the help of these 



two roots, we factorize the 3PN accurate expression for r 2 and obtain the following integral 
connecting t and s: 

f s - A + A 1 s + A 2 s 2 + A 3 s 3 + A 4 s 4 + A 5 s 5 
t-t = =^ ds. (9) 

The coefficients A t : % — 1..5 are some PN accurate functions of -E, /i and rj which follow 
directly from the expression for r 2 . The radial motion is uniquely parametrized by using the 
ansatz 

r = a r (l — e r cosw), (10) 

where a r and e r are some 3PN accurate semi-major axis and radial eccentricity which may 
be expressed in terms of s_ and s + as 

ls_ + s+ s_-s + 

a r = - , e r = ■ . (11) 

2 s_ s + s_ + s + 

With the aid of above relations and 3PN accurate s_ and s + we compute the 3PN accurate 
expressions for a r and e 2 in terms of E, h and 77. We now move on to obtain the 3PN 
extension to the 'Kepler equation'. For this purpose, we first compute the 3PN accurate 
expression for the radial period P associated with the radial motion. The expression for P 
directly follows from Eq. Q and reads 

p = 2 f s A + Ais + A 2 s 2 + A 3 s 3 + A4S 4 + A 5S 5 ^ ^ 
Js+ a/(s- - s)(s-s+) s 2 



eccentric anomaly u. We introduce an auxiliary variable v = 2 arctan 



Using Eq. 0, we express the mean anomaly / = n{t — t$) = -p-(t — to) in terms of the 

1/2 



tanf 

and with the aid of certain trigonometric relations involving v, displayed in Appendix [B 
we obtain the following temporary parametrization for the mean anomaly I which has the 
following general structure 

I = n(t — to) = u + kq sin u + -rriv — u) 

& 

+ — \ sin-u + -| sin2{; + sin3-u . (13) 
c 4 c 6 c 6 

The coefficients «j are some PN accurate functions of E, h and rj. Since these expressions 
are lengthy and are only required temporarily, we do not list them in the paper. 

Let us now turn our attention to the parametrization of the angular part, which is also 
required to obtain the final representation for the PN accurate 'Kepler equation'. Using the 
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3PN accurate expression for and the factorized 3PN accurate expression for s, we compute 
^ = 0/s and readily obtain the following relation between and s 

B + B 1 s + B 2 s 2 + B 3 s 3 + £? 4 g 4 + ff 5 s 5 ^ (M) 



V( s - - s )( s - s +) 

where the coefficients B{ : i = 1..5 are some PN functions of E, h and t] as presented in 
Appendix [X] Next step involves the computation of the advance of periastron during the 
radial period P. This quantity, denoted by $, is obtained by evaluating the integral on the 
right hand side of Eq. (fT4*|) between s + and s_, 

$ = 2 r b q + B ^ + B ^ 2 + ^ 3 + b 4 s 4 + £? 5 s 5 ^ 

7 S+ - s)(s-s+) 

To derive a temporary parametrization of the angular motion, using Eqs. (fT4"j) and (|15jh 
we first compute ^ (0 — 0o). Employing certain trigonometric relations, given in Appendix 
IB] we arrive at a temporary 3PN accurate parametrization for the angular motion 

— [0 — 0o J = « + -rSint) + - sin 2t> 
$ cr c 4 

+ —7 sin 35 H — ^sin4?)H — | sin 55, (16) 

c 4 c 6 c 6 

where A, are some PN accurate functions, expressible in terms of E, h and r\. 

To obtain the final parametrization for I and equations, we introduce v = 



2 arctan 



\ i/2 



tan^ 



where e</, differs from e r by yet to be determined PN correc- 



tions. It is straightforward to express v in terms of v, some trigonometric functions of v and 
the PN corrections relating and e r , as shown schematically in Eq. (jB2|) . We systematically 
introduce t> in Eq. (jTSj) such that the coefficient of sin w is zero up to the third PN order. This 
uniquely determines the PN corrections defining and consistently generalizes the imposed 
requirement that the expression for ^ (0 — 0o) at 1PN order has a 'Keplerian' structure. 
However, we find that at 2PN and 3PN orders, the v contribution will be supplemented by 
other trigonometric functions of v. With the introduction of v, the 3PN accurate expression 
for ^ (0 — O ), given by Eq. (jTfijl . becomes 

2tt , , 



$ - ~ M =V+ \ c* 



+ sin4tH — ^sin5w, (17) 
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where v = 2 arctan 



l+e 4 
1-ej 



1/2 



tan | 



. The explict 3PN accurate expressions for <E>, and 
the orbital functions will be displayed below. 

Let us now go back to Eq. (|13j). giving the parametrization to I in terms of u, e t and £(«). 
Using the 3PN accurate relation connecting v and v, we rewrite that equation and obtain 
the following 3PN extension to the 'Kepler equation': 



I = n (t — t ) 



9u 9et 
[ + J 

. fit . fet\ . . i>et . n . het . „ 
+ — - H — — sin u H — - sm 2 v H — sin 3 1> . 



(18) 



<r c u / c u c u 

The 3PN accurate expressions for n, e t and the orbital functions appearing in the above 
'Kepler equation' will be displayed below. 

Finally, we display, in its entirety, the third post-Newtonian accurate generalized quasi- 
Keplerian parametrization for a compact binary moving in an eccentric orbit in ADM-type 
coordinates 



r = a r (1 — e r cosw) , 
/ = n (t — to) = u — &t sin u + + ^f" J ( v ~ u ) 



, fit fet\ ■ . kt . h 6t . 
+ I — - H — — I sm v H — - sm I v H — — sm 3 u , 



(19a) 



(19b) 



2tt 



, /e<A • o , f 9H , 9e<f>\ . Q 



where v 
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+ — ^- sm4w H -r- smov . 
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L+e ^ ^ tan s 



(19c) 



l-e. 



The 3PN accurate expressions for the 
orbital elements e 2 , $, and e| and the post-Newtonian orbital functions 

9u, 9et, Ut, feu kt, het, fty, fe<f>, 9A<j>, 9e<t>, ie<p, and h 6( p, in terms of E, h and rj read 
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136 - 56r/ - 15r? 2 
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2x 



768c 6 (-2Eh 2 ) 2 

(-. 



31872 



-88404r/ + 2055r/7r 2 + 4176r/ 2 - 210r/ 3 + (^2256 

+ 10228r/ - 15r/7r 2 - 2406r/ 2 - 450r/ 3 ^ (-2Eh 2 ) + 677 (136 

+3477 + 31r/ 2 )(-2£/i 2 ) 2 J. (21b) 

These relations allow one to choose a specific eccentricity, while describing a PN accurate 
non-circular orbit. 

It is highly desirable that such a detailed parametrization and lengthy expressions for 
the PN accurate orbital elements and functions should be subjected to possible consistency 
checks. We have devised and performed the following consistency check for our computation. 
Note that the expressions for a T and e 2 were obtained from the PN accurate roots s_ and 
s+, representing the turning points of the radial motion. This prompted us to express the 
expressions for r 2 and <fi 2 , derived using the Hamiltonian equations of motion and given by 
Eqs. (JSJ), in terms of E, h, rj and (1 — e r cos(u)). We then compared the above expression for 
r 2 with the one derived using the parametrization, namely f 2 = (j- ^) 2 - This expression 
for ?" 2 , after some lengthy algebra, is found to be in total agreement with r 2 , computed 
using Hamiltonian equations of motion to the third post-Newtonian order. We note that 
the above computation fully checked both the structure and parameters of both the radial 
and temporal part of the generalized quasi-Keplerian representation. We performed a similar 
check on the angular part by noting that 2 = (^ 4^ ^) 2 . The expressions for 2 , computed 
using the above relation and via the Hamiltonian equations of motion, were also found to be 
in total agreement to 3PN order. These two computations provided us with powerful checks 
on our 3PN accurate generalized quasi-Keplerian parametrization. 

Finally, we observe that to the 2PN order, our results are in agreement with results 
available in 14, Q]. The above parametrization is also consistent with results given in 
[l^ . modulo typographical errors and omissions. In the next section, we obtain a similar 
parametrization in harmonic gauge. 
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IV. THE 3PN ACCURATE GENERALIZED QUASI-KEPLERIAN 
PARAMETRIZATION IN HARMONIC COORDINATES 



Incidentally, it was in harmonic coordinates that the quasi-Keplerian parametrization for 
compact binaries of arbitrary mass ratio moving in eccentric orbits was first realized jll|. 
The above computation employed 1PN accurate expressions for the conserved orbital energy 
and angular momentum in terms of the dynamical variables of the binary to derive 1PN 
accurate expressions, in harmonic coordinates, for r 2 and <j) 2 . However at 2PN order, the 



14, 



though the description of the 



parametrization was achieved only in ADM gauge 
binary dynamics, in these two gauges, begins to differ at 2PN order 

In this section, we will compute 3PN accurate generalized quasi-Keplerian parametriza- 
tion, in harmonic coordinates, for compact binaries in eccentric orbits. To begin our compu- 
tation, we will require 3PN accurate expressions for the conserved orbital energy and angular 
momentum in harmonic coordinates and in the center-of-mass frame. These quantities, writ- 
ten in terms of the dynamical variables of the binary, are required to derive, by iteration, 
3PN accurate expressions for r 2 and </> 2 in terms of E, h and r. The above mentioned 3PN 
accurate conserved orbital energy and angular momentum, in harmonic coordinates, are 



available in 



251 ] [For a comprehensive review of the post-Newtonian computations, espe- 



cially in harmonic coordinates, see |26( ]. However, we will not be able to employ directly 
these conserved quantities in our computation for r 2 and </> 2 for the following two reasons. 
First, at 3PN order, these expressions contain logarithmic terms involving the radial sepa- 
ration r. These logarithmic terms will not allow us to obtain 3PN accurate expressions for 
r 2 and <fi 2 as polynomials in s — 1/r, preventing the determination of 3PN accurate roots of 
r 2 and hence the parametrization in harmonic coordinates. However, in the near-zone of a 
gravitating system, the harmonic gauge conditions do not fix the coordinate system uniquely, 



allowing various harmonic coordinates to be employed at 3PN order 



27 



28 



2ai3o|. Further, 



it is possible to remove these logarithmic terms using a 3PN accurate coordinate transfor- 



mation that still respects the harmonic gauge condition |2j 



The problematic logarithmic 



N accurate expressions for the orbital energy and angular mo- 



terms appearing in the 3 
mentum, as available in |25[, were removed using the coordinate transformation given in 
Secondly, the expression for 3PN accurate orbital energy appearing in 25] involves an 
undetermined parameter, the so called regularization ambiguity, which was recently fixed by 
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different techniques (23 



llowed the compilation of following 
3PN accurate expressions for the orbital energy and angular momentum, in harmonic gauge 
and in the center-of-mass frame, reduced by fi as 



E = E + — Ei + — E 2 + E 3 , 

Cr C 4 C° 

J = |R x V| I J + Ji + -j J 2 + -« ^3 f 



(22a) 
(22b) 



cr c u 

where R and V are the relative separation and velocity vectors. The various contributions 
to E and J are given by 
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In the above expressions, r = |r| with r = R/ (GM), v 2 = v • v and r = — , where v = ^ 
and t is the coordinate time scaled by GM. The above expressions also match with those 
presented in Jgj]. 

Since the orbital motion is restricted to a plane, we express the components of r in polar 
coordinates as r = r (cos 4>, sin 4>) and obtain v = (r coscj) — r <t> sin^,?^ sine/) + r <j) cos<t>) 
allowing us to write v 2 
E and J, given by Eqs 



r 2 + r 2 2 . Using these inputs and 3PN accurate expressions for 



and 



we obtain, after lengthy iterations, 3PN accurate 



expressions for r 2 and 0, in terms of E, h = -^,f] and r. This leads to 3PN accurate 
expressions for r 2 and (d<p/ds), in harmonic coordinates, displayed in the appendix El We 
note that the expressions for r 2 and <f> differ from similar ones derived for ADM-type gauge 
at 2PN and 3PN orders. However, note that in both gauges the expressions for 7^ and <p 
are polynomials of degree seven. It is interesting to note that in the previous section the 
orbital dynamics was fully determined with the help of the Hamiltonian. However, here we 
employed neither the Hamiltonian nor the Lagrangian to determine the orbital dynamics. 
Though the energy is numerically equal to the Hamiltonian, we required both the energy 
and angular momentum as functions of particle positions and velocities to determine the 
binary dynamics. The analogy to energy as thermodynamical potential versus caloric and 
thermal equations of state is quite appropriate 3^|. 



17 



Following exactly the same procedure, detailed in the previous section, we derive the 3PN 
accurate orbital parametrization in harmonic gauge. The third post-Newtonian accurate 
generalized quasi-Keplerian parametrization, in harmonic coordinates, for compact binaries 
moving in eccentric orbits is given by 



r = a r (1 — e r cosu) 



I = n (t — to) = u — e t sin u + (J^- + ^-j (v — u) 
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where v — 2 arctan 



v + — - H 1- sin 2v 



c 4 c 6 



l^+^J sm3v 
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tan | 



The explicit 3PN accurate expressions for the orbital 
elements and functions of the generalized quasi-Keplerian parametrization, in harmonic 
coordinates, read 
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In harmonic coordinates too, there are PN accurate relations connecting the three eccen- 
tricities e r ,e t and e^. These relations read 
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We observe that the structure of the 3PN accurate parametrization is identical in both 
ADM-type and harmonic coordinates. This is not surprising as the expressions for f 2 and 
4> are polynomials of same degree in s = Mn these two gauges. We also observe that since 
the equations of motion, and hence the expressions for r 2 and 0, are the same in ADM 
and harmonic gauges to 1PN order, the orbital elements and the relations between the 
eccentricities to 1PN order are also identical in these two gauges. It is also possible to connect 
the radial eccentricities, and hence other eccentricities, associated with parametrizations in 
ADM-like and harmonic coordinates. This indicates that a circular orbit in ADM-type 
gauge also implies a circular orbit in harmonic gauge. The relation connecting e 2 in ADM- 
lik coordinates to that in harmonic coordinates reads 
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where e 2 |n and c 2 |a are the expressions e 2 in harmonic and ADM-type coordinates. 

The most striking result is that the expressions for n and $ in terms of E and h are 
identical to 3PN order in both ADM-type and harmonic coordinates. Indeed it was shown, 
using the Hamilton- Jacobi approach to describe the relative trajectory of the binary, that 
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the functional form of n and $ should be independent of the coordinate system used, if 
expressed in terms of gauge (coordinate) invariant quantities like E and h [l4| . This should 
be contrasted with the functional forms for other orbital elements like a r and e r , in terms 
of E and h, which depend on the coordinate system used. The explicit verification of the 
above prediction constitutes a powerful check on the algebra involved in the derivation of 
3PN accurate orbital representation in these two gauges. 

The gauge invariance of n and $ allows us to obtain following gauge invariant expressions 
for the third post-Newtonian accurate adimensional orbital energy and angular momentum 
for compact binaries of arbitrary mass ratio, moving in eccentric orbits. 
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The above equations generalize to eccentric 



orbits the gauge invariant expression connectin g th e orbital energy of a compact binary 
in circular orbit to its period at 3PN order [31I Isj] . Using 3PN accurate far- zone fluxes, 
which are not yet computed, and Eqs. (|28j). it should be possible to compute 3PN accurate 
expressions for dx/dt and dk'/dt in terms of x and k! . This would generalize the 3PN 
accurate gauge invariant expression for the rate of change of the orbital frequency of a 
compact binary in circular orbit, available in jj, [j^j]. The Eqs. (|28J) also indicate that it 
should be possible to characterize non-circular orbits in post-Newtonian relativity in an 
invariant manner. Therefore, this parametrization should be useful to analyze existing 
general relativistic simulations involving compact binaries in a post-Newtonian framework. 
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V. THE SUMMARY AND DISCUSSIONS 



In this paper, we have presented the third post-Newtonian accurate 'Keplerian' type 
parametrization for the motion of two non-spinning compact objects moving in an eccentric 
orbit. The above 3PN accurate parametrization, which is structurally quite similar to the 
2PN accurate generalized quasi-Keplerian representation, is given in both ADM-type and 
harmonic coordinates. The associated orbital elements and functions were explicitly com- 
puted in terms of the 3PN accurate conserved orbital energy, angular momentum and the 
finite mass ratio. We explicitly showed that to even these high post-Newtonian orders there 
are gauge invariant quantities to characterize an eccentric orbit in post-Newtonian relativity. 
We also performed, for the first time, some clever consistency checks to validate the lengthy 
algebraic and trigonometric manipulations involved in the derivation of the parametrization. 

There are quite a few possible applications for our parametrization and some of them 
are currently under investigation. The representation will be required to construct 'ready 
to use' search templates for inspiraling eccentric compact binaries, whose orbital dynamics 
is fully 3.5PN accurate, extending the currently available 2.5PN accurate ones, presented 
in P|. Recently, using a solution to 3PN accurate equations of motion for compact binaries 
moving in non-circular orbits, gauge dependent expressions for the associated conserved 
orbital energy and angular momentum were obtained which were used to analyze general 
relativistic simulations involving compact binaries [?], 0] • Using the arguments that there 
are gauge invariant expressions to 3PN order in our parametrization and that the orbital 
elements of the representation, when radiation reaction is included, are continuously evolving 
variables, it should be possible to construct a better post-Newtonian accurate 'diagnostic 
tool' to analyze the orbital configuration and motion of inspiralling compact binaries. 

Finally, we note that to obtain Oil/ c 11 ) corrections to the orbital averaged expressions 
for the far-zone fluxes associated with inspiralling eccentric binaries, 3PN corrections to 
the quadrapole approximation, our parametrization in harmonic coordinates will have to be 
heavily employed. The other crucial ingredients required for those computations, namely 
the 0(1/ c 11 ) corrections to the far- zone fluxes for compact binaries in general orbits, in 
harmonic coordinates, are expected to be available in the near future "\. 
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APPENDIX A: THE EQUATIONS GOVERNING THE 3PN ACCURATE RA- 
DIAL AND ANGULAR MOTION IN ADM-TYPE AND HARMONIC COORDI- 
NATES 



In this appendix, we display the 3PN accurate expressions for r 2 = j% and ^ = 
that are essential to obtain the generalized quasi-Keplerian parametrization in ADM-type 
and harmonic coordinates in terms of E, h, r] and s. As mentioned earlier, the expression 
for r 2 = (^f) 2 is required to obtain 3PN accurate turning points of the radial motion and 
hence required in the computations that determine the 3PN accurate radial parametrization 
and 'Kepler equation'. We first exhibit the 3PN accurate expression for r 2 in ADM-type 
coordinates. 
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To obtain the post-Newtonian accurate 'Kepler equation' and to parameterize the angular 
motion, we usually employ the PN accurate expressions for ^ factorized by the PN accurate 
roots s + and s_. 

Instead of displaying the original expression for g = 0/f , we exhibit below the expres- 
sion for factorized by the PN accurate roots s + and s_ which was actually employed in 
the computations. 
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In harmonic coordinates, the explicit expressions for r 2 and <fi are extracted from the 3PN 
accurate expressions for conserved orbital energy and angular momentum, in the center of 
mass frame, given by Eqs. (J22j) and (J23j) . To 3PN order, r 2 in harmonic coordinates reads 
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As in ADM-type gauge, the 3PN accurate factorized expression for reads 
d0 B + 5 lS + B 2 s 2 + B 3 s 3 + B 4 s 4 + 5 5 s 5 



ds v / (s--s)(s-s+) 
The explicit expressions for the £>j in harmonic coordinates are 
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We note that the expressions for t^ 2 and in these two gauges have the same structure 
though the coefficients differ at 2PN and 3PN orders. 



APPENDIX B: A SKETCH OF SOME COMPUTATIONAL DETAILS 



This appendix details the computation required to obtain the 3PN accurate 'Kepler 
equation' and the parametrization for the angular motion. As explained in Sec. IIH[ the 
temporary form for the 3PN accurate l(u) relation, given by Eq. (j!3j) . follows from Eqs. (jOJ) 
and (fT2*|) by employing the auxiliary variable v. This computation is not straightforward 
as the direct evaluation of the integral for (t — 1 ), given by Eq. (fT3|). gives a PN accurate 
expression in terms of E,h,r],S-,s + and s. To obtain the temporary form for the 3PN 
accurate 'Kepler equation', we multiply the above result with the 3PN accurate expression 
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for n, and use the trigonometric relations given below 

1 1 + e r cos v 

a r (1 — e r cos(u)) a(l — e 2 .) 



(Bla) 



M = arccos 2- — , (Bib) 

Vs_-s + (s_-s + )s/ 

( 2s s_ + s+\ 
v = arccos . (Blc) 

These relations are also employed heavily to obtain the temporary parametrization for the 
angular motion, given by Eq. (|16|). 

Let us turn our attention to the details of the computation, which gave the final 
parametrization for ^ ((f) — <p ). The starting points of the above calculation are Eq. (fTH|) 



1/2 



tan^ 



where 



and the introduction of PN accurate true anomaly v = 2 arctan , 
e^, differs from e r by yet to be determined PN corrections. It is easy to obtain the following 
3PN accurate expression for v in terms of v, which reads 

2 

y y 

v = v sin v + — r (sin 2v + 2 sin v) 

c 2 4c 4 

— y^~6 (3sin^ + 3sin2w + sin3f) (B2) 

where y is expressible in terms of those PN corrections which relate and e r . Using the 
above relation, we express ^g- (</> — </>o) > given by Eq. in terms of v and demand that there 
are no 'sinf ' terms to 3PN order. This requirement, as explained in Sec. lIII| is motivated by 
the desire that we want for the first post-Newtonian order 'Keplerian' like parametrization 
for the angular part. This procedure uniquely defines PN corrections connecting to e r 
and leads us to the final parametrization for the angular motion, given by Eq. (|17|) . 

Finally, we rewrite the temporary parametrization for the l(u) relation, given by Eq. (|T3*j) . 
in terms of v which leads us to the final expression for 3PN accurate 'Kepler equation' as 
given by Eq. (fTH|) . 
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FIG. 1: The geometrical interpretation of the eccentric and true anomalies, u and v, appearing 
in the Keplerian parametrization. The auxiliary circle of radius a circumscribes the orbital ellipse 
with semi-major axis a. The points O and F stand for the origin and the focus of the ellipse. 
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